5 research outputs found

    Perturbed, Entropy-Based Closure for Radiative Transfer

    Full text link
    We derive a hierarchy of closures based on perturbations of well-known entropy-based closures; we therefore refer to them as perturbed entropy-based models. Our derivation reveals final equations containing an additional convective and diffusive term which are added to the flux term of the standard closure. We present numerical simulations for the simplest member of the hierarchy, the perturbed M1 or PM1 model, in one spatial dimension. Simulations are performed using a Runge-Kutta discontinuous Galerkin method with special limiters that guarantee the realizability of the moment variables and the positivity of the material temperature. Improvements to the standard M1 model are observed in cases where unphysical shocks develop in the M1 model.Comment: 35 pages, 8 figure

    Models and numerical methods for time- and energy-dependent particle transport

    No full text
    Particles passing through a medium can be described by the Boltzmann transport equation. Therein, all physical interactions of particles with matter are given by cross sections. We compare different analytical models of cross sections for photons, electrons and protons to state-of-the-art databases. The large dimensionality of the transport equation and its integro-differential form make it analytically difficult and computationally costly to solve. In this work, we focus on the following approximative models to the linear Boltzmann equation: (i) the time-dependent simplified PN (SPN) equations, (ii) the M1 model derived from entropy-based closures and (iii) a new perturbed M1 model derived from a perturbative entropy closure. In particular, an asymptotic analysis for SPN equations is presented and confirmed by numerical computations in 2D. Moreover, we design an explicit Runge-Kutta discontinuous Galerkin (RKDG) method to the M1 model of radiative transfer in slab geometry and construct a scheme ensuring the realizability of the moment variables. Among other things, M1 numerical results are compared with an analytical solution in a Riemann problem and the Marshak wave problem is considered. Additionally, we rigorously derive a new hierarchy of kinetic moment models in the context of grey photon transport in one spatial dimension. For the perturbed M1 model, we present numerical results known as the two beam instability or the analytical benchmark due to Su and Olson and compare them to the standard M1 as well as transport solutions

    Models and numerical methods for time- and energy-dependent particle transport

    Get PDF
    Particles passing through a medium can be described by the Boltzmann transport equation. Therein, all physical interactions of particles with matter are given by cross sections. We compare different analytical models of cross sections for photons, electrons and protons to state-of-the-art databases. The large dimensionality of the transport equation and its integro-differential form make it analytically difficult and computationally costly to solve. In this work, we focus on the following approximative models to the linear Boltzmann equation: (i) the time-dependent simplified PN (SPN) equations, (ii) the M1 model derived from entropy-based closures and (iii) a new perturbed M1 model derived from a perturbative entropy closure. In particular, an asymptotic analysis for SPN equations is presented and confirmed by numerical computations in 2D. Moreover, we design an explicit Runge-Kutta discontinuous Galerkin (RKDG) method to the M1 model of radiative transfer in slab geometry and construct a scheme ensuring the realizability of the moment variables. Among other things, M1 numerical results are compared with an analytical solution in a Riemann problem and the Marshak wave problem is considered. Additionally, we rigorously derive a new hierarchy of kinetic moment models in the context of grey photon transport in one spatial dimension. For the perturbed M1 model, we present numerical results known as the two beam instability or the analytical benchmark due to Su and Olson and compare them to the standard M1 as well as transport solutions
    corecore